Inelastic shot noise characteristics of nanoscale junctions from first principles 
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We describe an implementation of ab-initio methodology to compute inelastic shot noise signals 
due to electron- vibration scattering in nanoscale junctions. The method is based on the framework 
of non-equilibrium Keldysh Green's functions with a description of electronic structure and nuclear 
vibrations from density functional theory. Our implementation is illustrated with simulations of 
electron transport in Au and Pt atomic point contacts. We show that the computed shot noise 
characteristics of the Au contacts can be understood in terms of a simple two-site tight-binding 
model representing the two apex atoms of the vibrating nano-j unction. We also show that the 
shot noise characteristics of Pt contacts exhibit more complex features associated with inelastic 
interchannel scattering. These inelastic noise features are shown to provide additional information 
about the electron-phonon coupling and the multichannel structure of Pt contacts than what is 
readily derived from the corresponding conductance characteristics. We finally analyze a set of Au 
atomic chains of different lengths and strain conditions and provide a quantitative comparison with 
the recent shot noise experiments reported by Kumar et al. [Phys. Rev. Lett. 108, 146602 (2012)]. 



I. INTRODUCTION 

The signatures of vibrational modes in the shot noise 
properties of nanoscale junctions have been the subject of 
active theoretical investigations [1-12]. Recently, it was 
shown in Refs. [6|-[8|, that under applying a bias volt- 
age eV larger than the typical phonon energy /icjo, the 
activation of phonon emission in a junction at low tem- 
peratures is responsible for a threshold behavior of the 
shot noise versus voltage characteristics. More specifi- 
cally, depending on the electronic transmission probabil- 
ity of the junction, the correction to the shot noise sig- 
nals induced by electron-phonon (e-ph) interactions was 
shown to exhibit jumps in the voltage derivative that 
are either positive or negative, as the result of a sub- 
tle interplay between one-electron tunneling events and 
correlated two-electrons processes [13]. This behavior of 
the inelastic shot noise signal was recently shown to be 
strongly dependent on fluctuations in the occupation of 
the locally excited vibrational mode. Under certain con- 
ditions, this phenomenon might lead to a strong feedback 
of the dynamics of the oscillator on the electronic noise 
properties p^, [ll] and the corresponding nonlinear effect 
in the shot noise could be of interest for the characteri- 
zation of heating effects at the nanoscale. 

In parallel to this theoretical activity, the first non- 
equilibrium shot noise measurements were recently per- 
formed on gold (Au) nano-j unctions that unravel clear 
signatures of the excitation of local vibrational modes 
P^. Those measurements confirm qualitatively the 
predictions of Ref. concerning the existence of a 

crossover from positive to negative correction to the Fano 
factor upon phonon excitation, when decreasing the elec- 
tron transmission factor r of the junction from unity. 
However, the experimental shot noise characteristics ex- 



hibit some unexplained features which seem to be out of 
the range of "single level, single vibrational mode" mod- 
els. For instance, the position of the crossover was shown 
to be shifted compared to the theory from r ~ 0.85 to 
r ?^0.95 [i3. 

The above mentioned experiments exemplify the rele- 
vance of exploring new methods that allow to compute 
quantitatively the inelastic shot noise signals from first 
principles. The aim of the present paper is thus two-fold: 
First, to document our implementation of such a frame- 
work into the Inelastica p^ - [T6| ab-initio code based 
on Siesta \yf\ and TranSiesta To this end, we 

adopt the results derived in Refs. 0, [20] based on non- 
equilibrium Keldysh Green's functions. Secondly, we ap- 
ply our implementation to discuss the inelastic shot noise 
signals in atomic point contacts of Au and Pt [21, 22] with 
all parameters extracted from atomistic calculations. De- 
pending on the material, a different number of conduc- 
tance channels are available for the electron transport 
[23l-[25| and — as a consequence — also the inelastic trans- 
port properties are qualitatively different [26] . In particu- 
lar, this allows us to highlight and analyze the additional 
complexities that arise from interchannel scattering un- 
der realistic conditions. We finally analyze calculations 
for a set of Au atomic chains of different lengths and 
strain conditions and compare the findings with the re- 
cent shot noise experiments [13]. 

The organization of the paper is the following. In 
Sec.ini we present the ab-initio methodology we have im- 
plemented in order to compute the correction to the shot 
noise induced by e-ph interactions. In Sec. lIIII we present 
results obtained for different Au and Pt atomic point con- 
tact, and calculations for Au atomic chains bridging the 
electrodes is investigated in Sec. IIVI Our conclusions are 
presented in Sec. [Vl 
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II. METHODOLOGY 

In this section we outline the methodology used to per- 
form first-principles calculations of shot noise character- 
istics. 



A. Model 

We consider the standard partitioning scheme in which 
an interacting device region D couples to two reservoirs 
of noninter acting electrons, namely left L and right R 
leads. This system is described by the following spin 
degenerate Hamiltonian 

H = Hd^Hl,r^Ht. (1) 

Here the device region in which the e-ph interactions 
are assumed to be strictly localized, is described by a 
Hamiltonian of the form 

Hn = Hi'^+HS!+H.-r>^, (2) 

HS! = T.f^AK (4) 

A 

^e-ph = E E ^i^-^'l^i(^A + h), (5) 

where and b]^ are the electron and phonon creation 

operators in the device space, respectively. Here ^i^^ 
is the single-particle Kohn Sham DFT Hamiltonian de- 
scribing electrons moving in a static arrangement of the 
atomic nuclei. In this Hamiltonian, electron-electron in- 
teractions are taken into account at the mean field level. 

is the phonon Hamiltonian of free uncoupled har- 
monic oscillators, and He-ph is the e-ph coupling within 
the harmonic approximation. 

The Hamiltonians describing the leads Hl^r and the 
tunnel couplings between leads and device region Hj^ are 
given by 

Hl,r= E J2^?i^Uc.,P (6) 

a—L,R ij 
a=L,R i,j 

where cj^ ■ is the electron creation operator in lead a = 
L, R. Each lead is considered to be in local equilib- 
rium such that the occupied states are characterized 
by a Fermi distribution with thermal energy ksT and 
chemical potential jUa- An applied voltage V is as- 
sumed to shift symmetrically the chemical potentials of 
the leads with respect to the Fermi level position at 
equilibrium Ep (determined self-consistently in Kohn- 
Sham DFT by filling the electronic states from below) 
as jiiL{R) = + (— )eV/2 and to leave the electrostatic 
potential of the device unchanged. 



B. Electronic structure methods 

All parameters in the above Hamiltonian are extracted 
from self-consistent calculations with Siesta [17] for the 
TranSiesta setup [3 according to the Inelastica 
scheme p^-[T6|. With {0^} denoting the full nonorthog- 
onal basis set of atomic orbit als, the Fermion opera- 
tors satisfy the anti-commutation relations {dl^dj} = 
{di^dj} = and {di^dj} = Sij, where Sij = 
is the overlap matrix (and similar for the lead Fermion 
operators). Using boldface notation throughout this 
manuscript for the electronic space, we let H^^^ repre- 
sent the matrix elements of the Kohn- Sham Hamiltonian 
in the device region, V" the coupling elements between 
device and lead a, S the overlap matrix, and M'^ the 
e-ph coupling matrix (obtained by finite differences [l5|) 
corresponding to a localized mode A with energy hujx- 

In the device region D the single-particle noninteract- 
ing retarded (advanced) Green's function g^^'^^(^), i.e., 
without e-ph interactions, takes the usual form 

gr(a)(^) = {{E ± io+)s - h(o) - 5:1^"^ (^) - j:'j^^\e)}-\ 

(8) 

where 

i:t\E) = {V"ygl^^\E)V" (9) 

is the retarded (advanced) self-energy due to lead a. Here 
g^J^^\E) represents the corresponding surface Green's 
function for the isolated lead and is calculated recursively 
[23. The retarded and advanced Green's functions are 
connected through the relation g^{E) = {g^{E)y . For 
convenience we also introduce here the level broadening 
due to lead a 

r„(i?) = i{s^(i?)-ss(i?)}. (10) 

Finally, as the theory of shot noise characteristics pre- 
sented in the following section is developed in an orthog- 
onal basis, we orthogonalize all relevant quantities ac- 
cording to standard Lowdin transformations [28], i.e., 

H(o) ^ = S-i/2h(0)S-i/2, (11) 

^ = S-i/^M^g-i/^^ (12) 

S ^ S = S-i/^SS-i/^ = 1, (13) 

g'-(<i) ^ gHa) = s-^/^g''''')S^/^, (14) 

r, ^ f „ = s-i/2r„s-i/2, (15) 

where the transformation matrices 

8^/2 = udiag(v^, . . . , v^)U-\ (16) 
= Udiag(l/v^, . . . , l/v^)U-i, (17) 

are determined through the eigenvalue problem SUi = 
SiUi for the overlap matrix. 
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C. Nonequilibrium Keldysh Green's functions 

We summarize in this section how to calculate conduc- 
tance and shot noise characteristics for nano-j unctions 
described by the electronic structure as outlined in 
Sec. IIIBI The concepts underlying this methodology are 
derived in Ref. [19] and were applied to single- level and 
single- mode models in Refs. [6|-l8|. Here we follow the 
specific generalization to the multi-level and multi-mode 
model formulated by Haupt et. al. [s*, ^205 . The two fun- 
damental approximations are (i) weak e-ph interactions 
and (ii) the so-called extended wide-band limit (EWBL). 

1. Extended wide-band limit (EWBL) 

As in previous work we adopt the EWBL [1, [H, [H, [20|, 
[29I . [30} to perform energy integrations analytically such 
that explicit results for the mean current and shot noise 
can be stated. The EWBL consists of approximating the 
noninteracting retarded and advanced Green's functions 
as well as the level broadenings with their values at the 
Fermi energy Ep^ i.e., 

g'-('^)(^)«g'^(»)(Ef) = (18) 

r«(^) «r„(^f) = r«. (19) 

Physically this is motivated by the fact that in many real 
systems the electronic spectral properties typically vary 
slowly on the scale of a few phonon energies and applied 
voltages p^, ll5|, |30| . In the case of atomic gold wires 
this approximation was successfully tested by one of us 
in Ref.[15| via a direct comparison to computationally 
more expensive calculations based on the self-consistent 
Born approximation. The physical reason is the strong 
hybridization of device states with the electrode states 
(life-time broadening on the eV scale) and the existence 
of only low-energy vibrational modes (on the meV scale). 
As this situation also applies to point contacts of Au and 
Pt, we expect the EWBL to be a very good approxima- 
tion for these systems too. 

2. Computing the current characteristics 

In absence of e-ph interactions (M^ = 0), the (bare) 
current Iq{V) is given by the standard Landauer-Biittiker 
formula [31]. Within EWBL the transmission function 
becomes energy independent so that the expression for 
the current- voltage characteristics is simply 

h[2e/h]{V) = Ti{T}eV, (20) 

where 

T = r^g^r^g^ (21) 

is the elastic (bare) transmission matrix of the junction. 
This expression is similar to the Fisher and Lee formula 
for the conductance 1321 . 



In presence of weak coupling to the vibrational sub- 
system (M'^ 7^ 0), the above expression for the mean 
electronic current has to be modified. At second or- 
der of perturbation theory in the e-ph coupling strength, 
the correction to the current 5I(y) (within EWBL) can 
be expressed in terms of products of microscopic factors 
(system dependent) by voltage-dependent universal func- 
tions (system independent) [15, 20]. We write the inelas- 
tic corrections to the current as 

5I[2e/h](y) = 6Iel{V) + (5/inel(^), (22) 

SU2e/h]{V) = E {(1 + 2n^)Tr{Tf } (23) 

A 

+2n^Tr{TW}}ey, 
6Iin42e/h]{V) = ^Tr{Tli)}5,(ey). (24) 

A 

with the microscopic factors given by traces over the fol- 
lowing quantities 

Tf = (g'-M^g^M^AH + H.c.) , (25) 
tI,'^ =rLg'^{M^AflM^ (26) 

-^(M^AMV^i^-H.c.)}g^ 

and the voltage dependent universal functions given by 

gx{eV) = eV^ l{U{eV - hux) - U{eV + hcux)}, (27) 
U{eV) = eFcoth(/3eF/2). (28) 

In the above equations, g^^ denotes the real part of g^, 
P = l/ksT the inverse temperature of the electrodes, 
and Aq, = g^rQ,g^ the partial spectral function corre- 
sponding to lead a. The total spectral function is given 
by A = Al + Ar. 

We note that the expressions Eqs. ([23]) - (|24|) are subject 
to the following three additional approximations: (i) We 
have ignored the asymmetric contributions to the con- 
ductance (with respect to voltage) which are derived in 
Ref. [20I, [33I, [33 • These contributions are logarithmically 
divergent in the zero-temperature limit for \eV\ = hux 
and signal the breakdown of second-order perturbation 
theory at the inelastic threshold. However, a resumma- 
tion scheme might renormalize and cure this problem 
jlO| . Out of the threshold region, the logarithmic cor- 
rections to the conductance are typically orders of mag- 
nitude smaller than the symmetric contribution. They 
also vanish in the limit of symmetric couplings to the 
left and right lead [15], thus justifying our assumption, 
(ii) We fix the phonon populations to the equilib- 
rium values as given by the Bose-Einstein distribution 
= l/[ex.p{f3huJx) — M (regime of equilibrated phonons). 
For atomic point contacts and chains this is a reasonable 
starting point as the vibrations in the nanoscale contact 
are damped to some extent by coupling to bulk phonons 
|l4 I15L l35l l36l] . Furthermore, the effect of phonon heat- 
ing on the shot noise characteristics is a delicate research 
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topic that is beyond the scope of the present study [lO|- 
[l^ . (iii) Equation (|22]) corresponds to taking into ac- 
count only the Fock (exchange) self-energy. Neglecting 
the Hartree term is indeed a good approximation in most 
cases of interest: within the EWBL this term provides 
a voltage independent renormalization of the molecular 
level positions in the regime of equilibrated phonons and 
thus displays no features at the phonon emission thresh- 
old (20j. 

The total correction to the current in Eq. (|22]) is the 
sum of two contributions. The first one 5Iqi{V) is due to 
elastic processes induced by e-ph interactions that renor- 
malize the bare transmission factor r to an effective one, 
e.g., f ~ r + Tr{T^^^} at zero temperature. The sec- 
ond contribution dIine\{V) originates from inelastic pro- 
cesses activated by phonon emission. The voltage de- 
pendence of this contribution is a universal function of 
voltage ^a(^) [see Eq. (|27|) ] that exhibits a threshold at 
|eF| = hux^ in the zero temperature limit. The sign and 
size of this threshold (jump in conductance) is controlled 
by the inelastic microscopic factors T^^^ given in Eq. (|26|) . 

Upon differentiation of Eq. (|22]) with respect to voltage 
one obtains the correction of the conductance dG{V) = 
dv{^I{V)). In the zero temperature limit, 6G{V) is dis- 
continuous at the inelastic threshold, due to the contri- 
bution of the inelastic term, see Eq. (|24|) . We thus define 
the corresponding jump in the inelastic correction to the 
conductance AGx = lim^^Q+ {6 G{huJx-\-r])—SG{huJx—r])} 
at the threshold voltage corresponding to mode A by 

AGA[2eV/i]=TV{TW}. (29) 

3. Computing the shot noise characteristics 

Shot noise characteristics in absence of coupling to lo- 
cal vibrational modes (M'^ = 0) was reviewed by Blanter 
and Biittiker in Ref. [37]. Within EWBL, the correlation 
function of the current operator evaluated at zero fre- 
quency, e.g. the (bare) shot noise characteristics 5*0 (^), 
is given by the simplified expression [s^ 

So[2eyh]{V) = ^Tt{T^} + Tr{T(l - T)}U{eV). (30) 

This expression is associated with the noise induced by 
thermal fluctuations in the electron occupation of the 
electrode Fermi seas and with fluctuations in the occu- 
pation of the coherent left- and right-moving scattering 
states. 

At second order of perturbation theory in the e-ph cou- 
pling strength, the expression for the finite temperature 
correction to the noise 5S{V) in the regime of equili- 
brated phonons is rather complicated, and derived in de- 
tail in Ref. [20j . While we have implemented these results 
and used them in the numerical part of this work, we 
here just state the simpler result for the zero-temperature 
limit 



5S[2eyh]{V) = SS.iiV) + SSinei{V), (31) 
^^eipeV/iKV') = ETr{(l-2T)Tf }|ey|, (32) 

A 

6Sinei[2eyh]{V) = ^TV{(1 - 2T)tW + Q^,} (33) 

A 

x{\eV\-nuJx)e{\eV\-fuJx), 

with 

Qa = -g'Tig'-{M^AKriAfiM^ (34) 
+M^AflrLg'^M^g'^rfl + H.c.}. 

Consistent with our assumptions for the inelastic correc- 
tions to the mean current f Sec. Ill we neglect also for 
the noise part the asymmetric terms leading to logarith- 
mic divergences as well as contributions from the Hartree 
e-ph self-energy. 

Analogous to the corrections to the current [Eqs. (|22]) - 
([28]) ] also the inelastic noise corrections [Eqs. (|3T]) - (|M|) in 
the zero-temperature limit] can be written as products of 
microscopic factors by universal voltage-dependent func- 
tions. The first term dSei{V) [Eq. (|32]) ] represents an 
elastic correction to the noise. The second term ^S'inei(^) 
[Eq. (|33]) ] is related to inelastic signatures of phonon ac- 
tivation in the shot noise. Its origin and interpretation is 
less intuitive than the corresponding expression for the 
mean current. The part proportional to Tr{(l— 2T)T^^'^} 
originates from a mean-field contribution to the shot 
noise while the other part proportional to Tr{QA} is 
related to vertex corrections [20]. A similar (although 
slightly different) decomposition in terms of one-electron 
(mean-field like) and two electron (vertex-like) processes 
was proposed in Ref. 

Instead of looking directly at the inelastic corrections 
in the shot noise it is convenient to analyze the voltage 
derivative of the shot noise SS{V) = dv{SS{V))^ i.e., 
the inelastic noise change. In the zero temperature limit, 
5Sx{V) is discontinuous at the inelastic threshold, due 
to the contribution of the inelastic term, see Eq. (|33]) . 
We thus define the corresponding jump in the inelastic 
correction to the shot-noise ASx = lim^^o+{^^(^^A + 
77) — SS{Jiwx — ^)} at the threshold voltage corresponding 
to mode A by 

ASx[2eyh] = Tr{(l - 2T)T^'^ + Qa}. (35) 

III. RESULTS FOR AU AND PT CONTACTS 

We first consider Au and Ft atomic point contacts, 
such as those shown Fig. [TJ as benchmark systems for 
our ab-initio calculations. Along the lines of Ref. [38j by 
one of us, we consider periodic supercells with a 4x4 rep- 
resentation of either Au(lOO) and Ft (100) surfaces sand- 
wiching two pyramids pointing toward each other. The 
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(a) 




Au contact 



FIG. 1. (Color online) Atomic point contacts of (a) Au atoms 
(shown in yellow-gold) and (b) Pt atoms (shown in blue-gray) 
considered in the first-principles transport calculations. The 
characteristic electrode separation L is measured between the 
second-topmost surface layers. The distance d characterizes 
the separation between the two apex atoms. 



characteristic electrode separation L is measured between 
the second- topmost surface layers, since the surface layers 
themselves are relaxed and hence deviate on the decimals 
from the bulk values. 

Our Siesta calculations use a single-^ plus polariza- 
tion (SZP) basis with a confining energy of 0.01 Ry 
[corresponding to the hd and 6(5, p) states of the free 
atoms], the generalized gradient approximation (GGA) 
for exchange-correlation, a cutoff energy of 200 Ry for the 
real-space grid integrations, and the F-point approxima- 
tion for the sampling of the three-dimensional Brillouin 
zone. The interaction between the valence electrons and 
the ionic cores is described by standard norm-conserving 
Troullier-Martins pseudopotentials generated from rela- 
tivist ic atomic calculations. As for the bulk Au and Pt 
crystals, we set the lattice constants to 4.18 A and 4.02 
A, respectively. 

For each electrode separation L we relax the surface 
atoms until the residual forces are smaller than 0.02 eV/A 
and proceed calculating vibrational modes and e-ph cou- 
plings by finite differences. For simplicity, we here only 
consider that the two apex atoms can vibrate, leaving us 
with six characteristic vibrational modes in the device. 
This assumption is only made to facilitate a fundamen- 
tal understanding of the inelastic signals. Finally, while 
electron transport in the supercell approach generally in- 
volves a samphng over -points, we approximate in the 
following all relevant quantities with their values at the 
F-point. 

We present in Fig. [2] the dependence of the electron 
transmission and vibrational frequencies, as a function 
of the electrode separation L. In both cases of Au and 
Pt junctions, the total transmission r decreases with 
L as observed in Fig. [2jb),(d) and covers the range 
from contact (ballistic limit) to the tunnel regime (low- 
transmission regime). The total transmission r = '^i 
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FIG. 2. (Golor online) (a) Vibrational frequencies and (b) 
electronic transmission at the Fermi energy as a function 
of electrode separation L for the Au atomic point contact. 

(c) ,(d) similar as (a),(b) but for the Pt atomic point con- 
tact. As indicated in the legends in panel (a) the longitudinal 
(transverse) eigenmodes are pictured in (a) and (c) with cir- 
cular (diamond) symbols. The out-of-phase (in-phase) modes 
are represented by filled (open) symbols. In panels (b) and 

(d) both the total transmission r = rtot (large red circles) as 
well as the transmission n , . . . , T4 for the four most conduct- 
ing transport eigenchannels (small symbols) are shown. Note 
that the legends shown in the Au panels apply also to the Pt 
panels. 



can be understood as a sum over eigenchannel trans- 
missions Ti for a set of (non mixing) electron scattering 
states. In Fig. [3] we have visualized the scattering states 
belonging to the four most transmitting channels (waves 
incoming from below) [39]. 

In the case of Au junctions the total transmission is 
essentially made up of a single channel, i.e., r ~ ri as 
seen in Fig. [2fb). This fact can be traced back to the 
single 5- valence of Au [111, [23l - [25| . The corresponding 
eigenchannel scattering state ipi is rotationally symmet- 
ric (cr- type) as seen in Fig.[3fa). The fact that essentially 
only one transmission channel contributes to the elastic 
current can also be appreciated by comparing the ampli- 
tude of the transmitted part of the different scattering 
states as it reflects the transmission probability. For the 
Au contact it is clear that only penetrates significantly 
the tunnel gap between the apex atoms. 
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FIG. 3. (Color online) Isosurface representations of the four most conducting eigenchannel scattering states V^^r) (incoming 
from below) for (a) a Au atomic point contact (L = 16.58 A) and (b) a Pt atomic point contact (L = 15.50 A). The states are 
ordered according to decreasing transmission. Due to the tunnel gap between the two sides the electron scattering states decay 
rapidly and the transmitted part of the wave on the other electrode is not visible. The blue and red colors represent the sign 
of the real part of the scattering states (our choice of phase makes the imaginary part negligible for visualization purposes). 
The isosurface plots reveal different rotational symmetry around an axis connecting the apex atoms. For both contacts one 
observes that the channels ipi and 1^4 are cr-type states (rotationally symmetric) while and are 7r-type states (with a 
nodal plane through the symmetry axis). 



For the Pt junctions on the other hand the total trans- 
mission has significant contributions from three eigen- 
channels in the contact regime, cf. Fig.[2fd). As revealed 
in Fig. [3l^b) the symmetry of the most transmitting chan- 
nel is of cr-type while the following two are of 7r-type with 
a nodal plane through the symmetry axis. This mul- 
tichannel nature reflects the partially filled sd valence 
shells of Pt [2i|,[23H25|. 

A. Au contacts 

1. 6G{V) and 6S{V) characteristics 

Using the methodology presented in Sec. Ill CI we pro- 
ceed by studying the inelastic effects in the transport 
through the considered Au atomic point contacts. Fig- 
ure m shows the curves obtained for the 6G{V) and 
6S{V) = dv{SS{V)) characteristics upon phonon exci- 
tation, for several electrode distances spanning the range 
from tunnel to contact. 

As shown in Fig. Hl^a), for each of the considered 



geometries the correction to the reduced conductance 
SG{V) exhibits a threshold-like character around \eV\ ~ 
huj^ y ^ 10 meV corresponding to the out-of-phase lon- 
gitudinal vibrational mode [38]. The signals from the 
other five modes are so small that they are hardly visi- 
ble. For the tunneling setups (r < 1/2) the activation of 
phonon emission processes above the inelastic threshold 
opens a new channel for conduction, thus increasing the 
conductance compared to its elastic background value. 
Contrary, in the contact regime (r > 1/2) the activation 
of inelastic scattering processes reduces the conductance 
(backscattering), i.e., results in a negative jump. 

Figure HJ^b) shows the corresponding SS{V) = 
dv{^S{V)) characteristics, which also exhibit a threshold 
response for voltages close to the phonon frequency of the 

^" mode. But in contrast to the conductance curves, 

the finite temperature effect is not only to smoothen the 
jump but also to produce some small downturn in the 
vicinity of the inelastic threshold [l, |20| . The sign of the 

jumps ASi y is consistent with the predictions based 

on "single- level, single- mode" models in Refs. [6|-[8|, i.e., 
that it is negative only in the region 0.15 < r < 0.85 and 
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FIG. 4. (Color online) Inelastic conductance and noise cor- 
rections for Au atomic point contacts with different electrode 
separations L in the regime of equilibrated phonons. (a) Con- 
ductance corrections 6G/t{V) induced by e-ph interactions as 
a function of voltage V. (b) Derivative of the shot noise with 
respect to voltage SS/r{V) induced by e-ph interactions. For 
each geometry six eigenmodes are considered as only the two 
apex atoms are vibrating, cf. Fig. [2ja). The calculations are 
performed at T = 4.2 K. 



positive elsewhere. The onset of a negative inelastic cor- 
rection to shot noise is not particularly intuitive. It was 
recently observed experimentally in shot noise measure- 
ments performed on Au nano-j unctions and explained 
in terms of correlated two-electron processes mediated 
by Pauli principle (Pauli blocking) and e-ph interactions 

A more direct way to appreciate these trends is shown 
in Fig. [5fa). Here the total jumps in the conductance 
AG = (blue circles) and derivative of shot noise 
versus voltage = AS'a (red diamonds) [over a 
voltage range that covers all possible phonon excitations] 
is shown as a function of the bare transmission r of each 
considered geometry. For comparison the correspond- 
ing jumps AGi y and A^^ y associated with just the 

most active ^" mode is shown with stars. All the 

computed data is consistent with sign changes in the in- 
elastic correction at r = 1/2 for the conductance and at 
r ~ {0.15,0.85} for the shot noise. 



2. Analytic model 

To develop an understanding for the calculated am- 
plitudes and sign changes for AG and AS shown in 
Fig. [5fa), we developed a simple two-site tight-binding 
model of the vibrating Au contact along the lines of 
Ref. ^] . This model is more appropriate to describe our 
physical problem than the single- level model analyzed in 
Refs. [sHl. 

As shown in the inset to Fig. [5fa) we represent each of 



FIG. 5. (Color online) Analysis of inelastic features for Au 
atomic point contacts with varying electrode separation d 
shown as a function of the total transmission factor r. (a) Ab- 
solute values of the total jump in the conductance AG (blue 
circles) and derivative of shot noise versus voltage AS (red 
diamonds) at zero temperature. Blue-dashed and red-dotted 
curves are analytic results for AG and AS, respectively, ob- 
tained within the corresponding two-site tight-binding model 
of the vibrating nano-junction (shown in inset), (b) The ra- 
tio AS/eAG (red circles) as a function of the transmission 
factor at zero temperature. The black-dashed curve represent 
the corresponding analytic result. Common to both panels: 
Stars correspond to the inelastic signals when retaining only 
the contribution of the longitudinal out-of-phase vibrational 

mode AG^ y and AS^ y. The parameters for the two-site 

tight-binding model are mo — 0.0167r, to — O.STSF, and 
eo — Ef (see text for details). 



the two apex atoms in the Au contact by a single orbital 
and write the noninteracting electronic Hamiltonian of 
this device region as 



H 



(0) 



eo t{d) 
t{d) eo 



(36) 



where eo is the onsite energy of each orbital (chosen to 
be equal to the Fermi energy Ep for simplicity) and t{d) 
is the hopping term that is modulated when varying the 
distance d between the electrodes. The hybridization of 
each orbital with its metallic electrode is described by 



F 






r 



(37) 



where F characterizes the coupling strength to each lead. 

For the dependence of t on distance d we adopt the 
simple relationship 



t{d) 



to 



l^e^d-do)/D^ 



(38) 



that interpolates between the tunneling regime {d:^ do) 
for which the hopping decreases exponentially with dis- 
tance and the contact regime {d ~ do) for which it de- 
creases linearly with distance. In Eq. (|38]) , the parameter 
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do is the typical distance where the contact is formed, to 
is an energy scale that provides the prefactor for the ex- 
ponential decay in the tunnel regime (taken to be larger 
than r/2) and D describes the size of the crossover region 
between the two regimes. 

Within the EWBL the elastic transmission for this two- 
site model is given by 



-(d) = ( 



rt{d) 



(r/2)2 + t(d)^ 



(39) 



As it is reasonable to consider T to be independent of 
d (and to be the largest energy scale of the model), we 
restrict the value of the hopping to be in the physically 
relevant branch < t{d) < T/2 < to, for which one 
sees that the correspondence between r{d) and t{d) is 
one to one, i.e., the transmission factor r(t) is a bijective 
function of the hopping spanning the range < r(t) < 1 
and simply decreases with d because of the dependence 
in Eq. (pj) . 

Finally, the two atoms are coupled to the out-of-phase 
longitudinal vibrational mode {< >) (vibrating at a fre- 
quency y) which modulates the interatomic distance 



d. The corresponding e-ph coupling matrix ^ can 
therefore be determined as the derivative of the electronic 
Hamiltonian with respect to i.e.. 



M^^[t{d)] mo 



t{d)[to - t{d)] 

r/2(to - r/2) [ 



1 

1 



(40) 



This expression for the e-ph coupling properly captures 
the physical behavior with namely ^ is propor- 

tional to t{d) in the tunnel regime (exponential depen- 
dence on hopping amplitude) and is almost constant in 
the contact regime (linear dependence on hopping am- 
plitude). The coupling matrix ^ depends on two 
parameters, namely mo which is the value of the e-ph 
coupling at unit transmission (obtained for t{d) = T/2) 
and the hopping energy scale to- The dependence of the 
coupling matrix with the distance d is being encoded into 
the one to one relation r[t{d)] and is thus not explicit in 
Eq. (|4Ql). 

In the zero temperature limit, the inelastic corrections 
to the mean current and shot noise for this two-site model 
can simply be expressed as 



SI[2e/h]{V) ^ 7(r; mo, to){2(l - r)eV + (1 - 2r){eV - hw^^)0{\eV\ - hw^^)^ 
dS[2eyh]{V) ^ 7(t; mo, to){2(l - r)(l - 2r)\eV\ + [l - 8r(l - r)] {\eV\ - huj^^)0{\eV\ 



4 



(41) 
(42) 



where 



, , ^ ( mo[to - t{T)] \ 

^("'"^°'^°)^H r/2(to-r/2) J 



(43) 



is an effective e-ph coupling constant that depends on all 
the parameters describing the electronic and vibrational 
structure of the junction, namely r, mo and to- 

The comparison of our ab-initio results with Eqs. (|4T]) - 
(|42]) is shown in Fig. [5fa). The blue-dashed and red- 
dotted lines correspond to the results for the analyt- 
ical jumps AG and A^, respectively. We modulated 
the transmission factor r by decreasing the value of the 
hopping term from t{d) = r/2 (in the contact regime 
r = 1) to t{d) ~ (in the tunnel regime r ~ 0). 
We found that a reasonable (although not perfect) fit 
to the ab-initio data points could be achieved by fix- 
ing the two independent parameters mo = O.OIGTF and 
to = 0.875r. It is interesting to notice that the simple 
analytical model of Refs. [6-8] fails to reproduce both 
the shape and amplitude of the curves in the full range 
of transmissions r G [0, 1] (not shown here), mainly be- 
cause the e-ph coupling strength changes when varying 
the distance between the electrodes in a way qualitatively 
provided by Eq. (|4Q|) . We also note that the analytic re- 
sults in Fig. [5fa) display clear asymmetries with respect 
to r = 1/2. This is because the situations r ^ (tunnel 
limit) and r ^ 1.0 (ballistic limit) are physically very dif- 



ferent (this asymmetry is also present in the single-level 
models of Refs. [6-Sj)- Moreover, the inelastic correc- 
tions, as given by Eqs. (|^T]) - (|^2]) - (|^3]) . have nontrivial 
dependences on r. For the chosen model parameters we 
find extrema in AG around r ~ {0.22,0.94} and in /S.S 
around r ^ {0.07,0.55,0.99}. 

A remarkable feature of the two-site tight-binding 
model is that the ratio of /^S[2e^/h] to eAG[2eV^] is 
a universal function of r and is independent of the effec- 
tive e-ph coupling strength 7(r;mo,to), ie.. 



/^S[2e^/h] 1 - 8r(l - r) 



eAG[2e7/i] 



1 -2r 



(44) 



This result is also found for the single-level model of 
Refs. as a common prefactor containing the details 
of the electronic structure cancels out. Figure [5fb) shows 
that our ab-initio data follows quantitatively the analytic 
results for the ratio AS/eAG (black-dashed curve). The 
agreement with the analytical result is even better when 
considering only the contribution from the out-of-phase 
longitudinal vibrational mode to the inelastic signals, i.e., 
AG^ — ^ and AS^ — ^ (red stars). 
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B. Pt-contacts 



1. 5G{y) and 5S{y) characteristics 
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FIG. 6. (Color online) Inelastic conductance and noise cor- 
rections for Pt atomic point contacts with different electrode 
separations L in the regime of equilibrated phonons. (a) Con- 
ductance corrections dG/riV) (in the contact regime) induced 
by e-ph interactions as a function of voltage V . (b) Derivative 
of the shot noise with respect to voltage 5S /t{V) (in the con- 
tact regime) induced by e-ph interactions, (c)-(d) as (a)-(b) 
but for the geometries corresponding to the tunneling regime. 
For each geometry six eigenmodes are considered as only the 
two apex atoms are vibrating, cf. Fig. [2jc). The calculations 
are performed at T = 4.2 K. 

We present in Fig. [6] the corresponding results for Pt 
atomic point contacts as was given in Fig. S] for Au con- 
tacts. In contrast to the Au case the electronic transport 
properties of Pt contacts can no longer be understood 
in terms of a single conducting eigenchannel. In fact, as 
seen from Fig.[2fd) the contact regime is characterized by 
three almost open channels, i.e., one a- type (labelled ) 
and two 7r-type (labelled V^J 3) visualized in Fig. [3fb). 
The fourth channel ^/^J is included as it turns out that 
scattering into such closed channels are important to un- 



derstand the inelastic transport characteristics. As for 
the Au contact the transmission factor decreases with 
L and drops suddenly at the point where the contact 
(chemical bond) breaks (r ~ 2.14 for the critical geom- 
etry L = 14.80 A). Beyond this point the transmission 
drops exponentially with d signalling the tunnel regime 
[see Fig.EJd)]. 

The case of the contact regime (r ~ 3.0) is shown 
in Fig. [6l^a)-(b). The inelastic features in the SG{V) 
and dS{V) characteristics reveal several steps associ- 
ated to the excitation of transverse and longitudinal vi- 
brational modes. In this transport regime, the jumps 
AG^^"^ {AS^^"^) are not necessarily negative (positive) — 
i.e., dominated by inelastic backscattering processes — as 
expected for a single-channel system close to the ballistic 
limit. As shown in Fig. [6fa)-(b), the sign of the jumps 
can be the other way around (see also the discussion in 
Sec. nil B 2p . When the transmission factor decreases, 
one enters into the tunnel regime (the junction breaks). 
As shown in Fig. [6l^c)-(d), two main inelastic signals are 
seen and the jumps AG(^) and AS'(^) are always positive 
in the case of low transmissions. 



2. Mode by mode analysis 

Due to the multichannel nature of the Pt contacts, 
the 6G{V) and SS{V) characteristics cannot be described 
within the framework of the simple analytical model used 
for Au contacts in Sec. IIII A 21 Instead, we can gain an 
understanding for the characteristics by analyzing the 
contribution from each vibrational mode to the total in- 
elastic signals dG{V) and SS{V) for the Pt contacts. For 
simplicity, we restrict our analysis to the L = 14.50 A ge- 
ometry representative for the contact regime [Fig. [7ja)- 
(b)] and to the L = 15.50 A geometry representative for 
the tunnel regime [Fig. [71^c)-(d)]. These two geometries 
are shown respectively in Fig. [6] with plain red and blue 
lines. 

Table H] lists the vibrational modes and energies hwx, 
the corresponding inelastic corrections in conductance 
AGx and shot noise ASx, and a decomposition of the 
underlying scattering processes among the eigenchannels. 
The idea is that inelastic scattering can be understood in 
terms of Fermi's golden rule where scattering occurs from 
occupied eigenchannel scattering states i/j^ (for channel 
i) originating in the left electrode into empty scatter- 
ing states ipj^ (for channel j) originating in the right 
electrode (or vice versa, depending on the bias polar- 
ity) [4dl, [4l|. The total scattering rate, proportional to 
j KV^f l^'^l'^j^)!^' therefore be decomposed into 
intrachannel {i — j) and interchannel {i ^ j) compo- 
nents. 

We begin our analysis with the results for the contact 
geometry {L = 14.50 A) as shown in Fig. [7fa)-(b). The 
conductance 5G{V) curve is dominated by the inelastic 
contribution from the longitudinal, out-of-phase vibra- 
tional mode (labelled < >) with an energy quantum of 
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FIG. 7. (Color online) Mode by mode analysis of inelastic 
features for two characteristic Pt atomic point contacts in the 
contact and tunneling regimes, respectively, (a) Correction 
to the conductance 5G/t and (b) correction to the derivative 
of the shot noise with respect to voltage 5S/t from the six 
characteristic vibrational modes as a function of voltage eV 
for L — 14.50 A (contact regime), (c)-(d) as (a)-(b) but for 
L — 15.50 A (tunnel regime). The calculations are performed 
at T — 4.2 K. The mode character is shown for each curve 
with two arrows similar to Fig.[2{c). The total characteristics 
5G/t{V) and 5S /riV) (sum over all vibrational modes) are 
shown as colored plain lines. 



hjOi y = 16.2 meV. As reported in Tab. HI AG^ ^ is 

found to account for 116.7 % of the overall conductance 
correction. The table also reports that for this partic- 
ular mode intrachannel inelastic transitions are clearly 
dominant. In fact, 88 % of the total scattering corre- 
sponds to intrachannel scattering involving the 10 most 
transmitting eigenchannels. Furthermore, transitions be- 
tween the cr-type and 7r-type states shown in Fig. [3] are 
symmetry forbidden (0 % scattering) as expected for a 
longitudinal mode that creates a rotationally symmetric 
deformation potential along the transport axis. The e-ph 
coupling matrix ^ is therefore essentially diagonal in 
the eigenchannel basis. Finally, since scattering from this 



Mode A 


i > x2 tt x2 


fiuj\ [meV] 




16.2 


11.6 


9.8 6.7 






116.7 


38.4 


-38.8 -16.3 


A5'a[%A5'] 




72.0 


43.0 


51.1 -66.1 


intra^: V^f ^ V^f [i 


< 10) [%] 


88 


61 





inter^: V^f ^ {i 


^J) [%] 


11 


34 


97 92 


inter^: V^f,4 ^ V^J,3 


[%] 








70 66 



(b) L = 15.50 A - tunnel regime 



Mode A 


i > x2 tt x2 


hjjx [meV] 




15.7 


16.3 


12.0 11.8 


AG'a[%AG'] 




65.4 


-0.1 


26.2 8.5 


A5'a[%A5'] 




61.1 


-0.1 


29.4 9.6 


intra^: ^ V^f {i 


< 10) [%] 


99 


36 





inter^: ^ {i 


^3) [%] 





46 


99 97 


inter^: V^f,4 ^ V^?,3 


[%] 








65 27 



^ considering the 10 most transmitting eigenchannels 
^ considering only scattering among eigenchannels 1-4 

TABLE I. Mode by mode analysis of the vibrational frequen- 
cies hjj\ and jumps in the inelastic signals {/\G^ and A5''^) 
for two Pt atomic point contacts. The characteristic elec- 
trode separation is (a) L — 14.50 A (contact regime) and (b) 
L — 15.50 A (tunnel regime). Relative contributions from 
inter- and intrachannel scattering processes to the total scat- 
tering rate are given in percent. 



mode is essentially intrachannel involving only the three 
most transmitting eigenchannels, its effect can be ratio- 
nalized in terms of a simple summation over three inde- 
pendent single-channel models where the resulting 

jump AG^ y is expected to be negative and AS*^ ^ to 

be positive because the condition > 0.85 is satisfied for 
each of the first three eigenchannels, cf. Fig.[2fd). Indeed 
this is consistent with the numerics in Fig. [71^a)-(b). 

As for the remaining vibrational modes the contribu- 
tions to 5G{y) are rather small compared to the out-of- 

phase longitudinal mode {i >) as quantified in Tab. [J 

[42[. However, their contributions in the 5S{V) charac- 
teristics are — interestingly — much more pronounced. As 
shown in Fig.[7Jb) we can identify several inelastic signals 
in the shot noise of comparable order of magnitude that 
were not clearly visible in the conductance. An interest- 
ing case is the negative jump /S.S^^ due to the excitation 
of the in-phase transverse vibrational modes (tt? doubly 
degenerate). Due to symmetry this mode does not al- 
low for intrachannel scattering (0 %) and its effect can 
therefor not be rationalized in terms of single-channel 
models In fact, this highlights that the e-ph cou- 

pling matrix M^^ is essentially off-diagonal in the eigen- 
channel basis. Indeed the jump AS*^^ is negative despite 
that Ti > 0.85 for the three most transmitting eigen- 
channels, contrary to the understanding derived from a 
single-channel picture. The clear signature in the shot 
noise from the tt mode is therefore a prominent demon- 
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stration that information about the e-ph couphng may 
be extracted from noise measurements despite that the 
mode is essentiahy passive in the corresponding conduc- 
tance characteristics. 

We complete our analysis by considering also a repre- 
sentative case in the tunneling regime (L = 15.50 A) 
shown in Fig. Wic)-{d). The SG{V) and SS{V) char- 
acteristics both exhibit positive jumps of comparable 
heights. The dominant inelastic features are located 

at voltages corresponding to hw^ y = 15.7 meV and 

hhj^l = 12.0 meV that are associated with the excitation 

of the out-of-phase longitudinal mode {< >) and of the 

two degenerate out-of-phase transverse modes (ti x2), 

respectively. According to Tab. |I] the < > mode again 

gives rise only to intrachannel transitions (99% of the 
total scattering involves intrachannel scattering for the 
10 most transmitting eigenchannels). The e-ph coupling 
matrix ^ is therefore almost diagonal in the eigen- 
channel basis and the physics can again be understood in 
terms of a summation over independent channels. The 

fact that the jumps AG^ y and A^^ y are both positive 

is consistent with the picture from single-level models as 
the condition < 0.15 is satisfied for all eigenchannels, 
cf. Fig. [Sid). 

Considering the out-of-phase transverse modes (t| x2) 
we again find that the main inelastic processes are inter- 
channel transitions (see Table H]) and therefore that the 
e-ph coupling matrix M^^^^ is off-diagonal in the eigen- 
channel basis. In essence, transitions occur between the 
cr-type and 7r-type states shown in Fig. [3] as expected for 
a transversal mode due to symmetry. Again, single-level 
models are inadequate to describe the physics and there- 
fore do not predict the sign of the jumps AG^^ and AS^^. 
According to our numerics we find that both jumps are 
positive. 



IV. RESULTS FOR AU ATOMIC CHAINS 

Au atomic chains constitute benchmark systems for 
the exploration of inelastic signatures in the transport 
characteristics [15, 29, 30, 35, 43] and most recently the 
effect of phonon emission in the electronic shot noise was 
reported for the first time [13]. To complete our study 
we present in this section realistic calculations of inelas- 
tic signals in the SG{V) and SS{V) characteristics for 
a set of Au atomic chain geometries, shown in Fig. [8l 
corresponding to different lengths and strain conditions. 

Our calculations were carried out in the same way as 
described in SecHIIlfor the Au and Pt atomic contacts, 
except that now we allow all atoms bridging the Au(lll) 
surfaces to vibrate, i.e., the dynamical region consists of 
15 vibrating atoms. Again, for each electrode separation 
L we relax the dynamical atoms and the topmost Au(lll) 
layer until the residual forces are smaller than 0.02 eV/A. 
All the structures considered here are stable in the sense 
that hux > for all vibrational modes (no imaginary 
frequency modes exist in the calculations). 



In order for us to explore the location of a possible sign 
change in SS{V) we intended to generate very different 
geometries with transmission factors spanning as large 
an interval as possible. For all the considered structures, 
shown in Fig. [51 we obtained the range r = 0.911 — 1.02 
which is significantly more narrow than reported in the 
experiments of Ref. [13j . Within our treatment it appears 
difficult to construct geometries where the transmission 
deviates significantly from unity (unless the chain is bro- 
ken). 

The generic behavior for the Au atomic chains — such 
as the plain red curve for the geometry L = 27.00 A 
in Fig. [9l^a) (for which r = 0.991) — is a main inelas- 
tic threshold in the conductance curve dG{V) located at 
voltages \eV\ ^ 10 meV. The corresponding negative 
jump in conductance is due to inelastic backscattering 
events associated with the activation of the out-of-phase 
longitudinal (alternating bond length) vibrational mode 
of the atomic chain (the zone-boundary phonon mode 
with wave length q ^ 2/ci?) [lH, |43|. This active mode 
is sketched qualitatively with arrows in Fig. [51 The size 
of the jump is of order 1%, which is significantly larger 
than the corresponding signals in the case of Au contacts, 
cf. Fig. [ll^a). This is consistent with tight-binding models 
[sol, 13 which have shown that the inelastic conductance 
correction is roughly proportional to the number of vi- 
brating atoms in the chain. When considering the noise 
properties one observes a positive jump of similar mag- 
nitude in the SS{V) characteristics as seen in Fig. [9l^b). 

In order to compare the output of our numerical calcu- 
lations to the experimental results of Kumar et al. [1^ , 
we represented in Fig. Ws) the total shot noise charac- 
teristics S{y) = Sq{V)^5S{V) expressed in the reduced 
scale Y{V) = {S{V)-S{V = 0))/S{V = 0), as a function 
of the reduced parameter X{V) = l3eV/2coth{l3eV/2). 
As in Ref. [13], the Y{X) curves are well approximated 
by piecewise linear functions. When eV is below the 
main inelastic threshold, the slope of the Y{X) curve 
gives the Fano factor F_. Upon phonon excitation, the 
slope changes to thus defining a modified Fano factor 
(due to the inelastic correction to shot noise SSinei{V) 
[Eq. (|33l) ]) and a relative jump of Fano factor AF/F = 
(F+ - F_)/F_ (see Ref. [Ts]). We represent on Fig.[9i;d) 
the dependence of the relative jump of Fano factor AF/F 
(red points) as a function of the renormalized transmis- 
sion factor f ~ r + XIa '^^{'^a'^} (which coincides with 
the zero bias conductance). Each point is associated to a 
given geometry shown in Fig. [5] and is computed by nu- 
merical differentiation of the corresponding Y(X) curve. 
The experimental points from Kumar et al. are 
shown in the range r = 0.9 — l.Oas blue crosses. 

We first remark that for transmission factors r > 0.95 
the computed points agree quantitatively with the ex- 
perimental ones. This exemplify the predictive power of 
ab initio methods for generating quantitative results for 
transport calculations without any adjustable parameter. 
Although our results are fully consistent with the experi- 
mental data from Ref. [isl] in the regime r > 0.95, we are 
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FIG. 8. (Color online) Stable Au atomic chain geometries considered in the first-principles transport calculations. The 
characteristic electrode separation L, measured between the second-topmost surface layers, is indicated over each geometry. 
The active out-of-phase longitudinal (alternating bond length) vibrational mode is sketched qualitatively with arrows in one 



unable to identify a sign change in the correction to the 
shot noise, which was reported experimentally to occur 
around r 0.95 instead of the expected theoretical value 
r ~ 0.85 (see Fig. [9Kd) blue crosses). 

In order to try to extrapolate the position of the 
crossover as obtained in the output of our calculations, 
we summarize in Fig. [TOT a) the dependence of the size of 
the total jumps AG and as a function of the transmis- 
sion factor r of the Au atomic chains. This plot clearly 
shows that statistically the jumps in conductance AG 
(blue points) and shot noise A^ (red points) have oppo- 
site sign and very similar magnitude. The more detailed 
correlation between AG and AS* can be characterized by 
considering the ratio AS/eAG as shown in Fig. [TOlf b). 
Here it is observed that the absolute value of the ratio 
tends to increase with r. However, the computed ratios 
(red points) do not fall exactly on the prediction of an- 
alytical models (black dashed curve) such as the result 
in Eq. (jHj). Rather it appears as if the analytics is an 
upper bound to the computed ratios [45] . The data that 
agrees the best with the analytic ratio corresponds to the 



cases of even-numbered, linear atomic chains L = 27.00 
A and L = 27.20 A. 

Our data in Fig. [TOTb) actually suggests that the 
crossover could occur at a lower value than the r ^ 0.85 
value from analytic models. An interesting future exten- 
sion of our results would be to try to include disorder- 
induced conductance fluctuations in the treatment, as 
this effect played an important role in the interpreta- 
tion of the experiments in Ref. [13j . Such a development 
with a semi-empirical model of disorder could be achieved 
through a multi-scale methodology [46^, but at the cost 
of giving up the parameter-free scheme presented here 
for simulating the inelastic transport characteristics of 
nanodevices. 



V. CONCLUSIONS 

We have implem ented i nto the Inelastica p^-flGj 
code a methodology [19, 20] that allows to compute quan- 
titatively inelastic shot noise signals from first princi- 
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FIG. 9. (Color online) Inelastic conductance and noise char- 
acteristics for Au atomic chains shown in Fig. [8] in the regime 
of equilibrated phonons. (a) SG/r{V) and (b) SS/r{V) char- 
acteristics computed at T = 4.2 K. (c) Total shot noise 
S{V) = So{V) + SS{V) expressed in the reduced scale Y{V) = 
{S{V) - S{V = 0))/S{V = 0) as a function of the reduced pa- 
rameter X{V) = ^eV/2coth{l3eV/2) computed at T = 4.2 K. 
The dashed curves are the extrapolation of the low voltage re- 
duced noise characteristics Y-{X) = — 1), with F- the 
Fano factor (see the text for details), (d) The relative jump of 
Fano factor AF/F = [F^ — F-)/F- (red points) as a function 
of the renormalized transmission factor f ~ TrjT^^'*}, 
as extracted numerically from the Y[X) characteristics. For 
comparison, the blue crosses are the experimental results of 
Kumar et al. [13] . 



pies. The method was illustrated for Au and Pt atomic 
point contacts as well as Au atomic chains described 
at the atomistic level. We showed that the Au con- 
tact constitutes a benchmark system that can be un- 
derstood in terms of the "single level, single vibrational 
mode" models pJsj. To be quantitative we rationalized 
computed inelastic corrections to the conductance and 
shot noise characteristics in terms of a "two-site" tight- 
binding model for the vibrating Au contacts. However, in 



the case of Pt contacts the multichannel nature compli- 
cates the physics and simple analytic models were shown 
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FIG. 10. (Color online) Inelastic conductance and noise cor- 
rections for Au atomic chains shown in Fig. [8] in the regime 
of equilibrated phonons. (a) Total jumps AG (blue points) 
and AS (red points) at zero temperature as a function of the 
transmission factor, (b) The ratio AS/eAG (red points) as a 
function of the transmission factor. The black dashed curve is 
the prediction of analytical models for the ratio, cf. Eq. ([H]) . 



to be inadequate. Interestingly, we found that for Pt con- 
tact the coupling to transverse vibrational modes gives 
rise to features in the inelastic shot noise signals (origi- 
nating in interchannel scattering processes) that are not 
readily visible in the conductance. This opens up an al- 
ternative approach to characterize the e-ph couplings and 
the underlying multichannel electronic transport prob- 
lem. 

We also analyzed Au atomic chain configurations in 
order to compare directly with the experimental results 
of Ref. [13] for the inelastic effects in the shot noise. 
While the simulated shot noise behavior was found to 
agree quantitatively with the experiment in the ballistic 
regime (0.95 < r ^ 1) we were unable to identify the sign 
change in the shot noise correction that was reported in 
Ref. pj]. 

The techniques for characterizing inelastic effects in 
shot noise appears as a promising way to gain deeper 
insight into transport properties of nanoscale devices. It 
is our hope that the methodology presented in this paper 
will stimulate further research in this direction. 
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